Improving HLA typing imputation accuracy and eplet identification with local next-generation sequencing training data.

Assessing donor/recipient HLA compatibility at the eplet level requires second field DNA typings but these are not always available. These can be estimated from lower-resolution data either manually or with computational tools currently relying, at best, on data containing typing ambiguities. We gathered NGS typing data from 61,393 individuals in 17 French laboratories, for loci A, B, and C (100% of typings), DRB1 and DQB1 (95.5%), DQA1 (39.6%), DRB3/4/5, DPB1, and DPA1 (10.5%). We developed HaploSFHI, a modified iterative maximum likelihood algorithm, to impute second field HLA typings from low- or intermediate-resolution ones. Compared with the reference tools HaploStats, HLA-EMMA, and HLA-Upgrade, HaploSFHI provided more accurate predictions across all loci on two French test sets and four European-independent test sets. Only HaploSFHI could impute DQA1, and solely HaploSFHI and HaploStats provided DRB3/4/5 imputations. The improved performance of HaploSFHI was due to our local and nonambiguous data. We provided explanations for the most common imputation errors and pinpointed the variability of a low number of low-resolution haplotypes. We thus provided guidance to select individuals for whom sequencing would optimize incompatibility assessment and cost-effectiveness of HLA typing, considering not only well-imputed second field typing(s) but also well-imputed eplets.

Imaging-guided prognostic score-based approach to assess the benefits of combotherapy versus monotherapy with immune checkpoint inhibitors in metastatic MSI-H colorectal cancer patients.

BACKGROUND: This retrospective study determined survival responses to immune checkpoint inhibitors (ICIs), comparing mono- (mono) and combo-immunotherapy (combo) in patients with microsatellite instability-high (MSI-H) metastatic colorectal cancer (mCRC) by analyzing quantitative imaging data and clinical factors. METHODS: One hundred fifty patients were included from two centers and divided into training (n = 105) and validation (n = 45) cohorts. Radiologists manually annotated chest-abdomen-pelvis computed tomography and calculated tumor burden. Progression-free survival (PFS) was assessed, and variables were selected through Recursive Feature Elimination. Cutoff values were determined using maximally selected rank statistics to binarize features, forming a risk score with hazard ratio-derived weights. RESULTS: In total, 2258 lesions were annotated with excellent reproducibility. Key variables in the training cohort included: total tumor volume (cutoff: 73 cm3), lesion count (cutoff: 20), age (cutoff: 60) and the presence of peritoneal carcinomatosis. Their respective weights were 1.13, 0.96, 0.91, and 0.38, resulting in a risk score cutoff of 1.36. Low-score patients showed similar overall survival and PFS regardless of treatment, while those with a high-score had significantly worse survivals with mono vs combo (P = 0.004 and P = 0.0001). In the validation set, low-score patients exhibited no significant difference in overall survival and PFS with mono or combo. However, patients with a high-score had worse PFS with mono (P = 0.046). CONCLUSIONS: A score based on total tumor volume, lesion count, the presence of peritoneal carcinomatosis, and age can guide MSI-H mCRC treatment decisions, allowing oncologists to identify suitable candidates for mono and combo ICI therapies.

A novel batch-effect correction method for scRNA-seq data based on Adversarial Information Factorization.

Single-cell RNA sequencing (scRNA-seq) technology produces an unprecedented resolution at the level of a unique cell, raising great hopes in medicine. Nevertheless, scRNA-seq data suffer from high variations due to the experimental conditions, called batch effects, preventing any aggregated downstream analysis. Adversarial Information Factorization provides a robust batch-effect correction method that does not rely on prior knowledge of the cell types nor a specific normalization strategy while being adapted to any downstream analysis task. It compares to and even outperforms state-of-the-art methods in several scenarios: low signal-to-noise ratio, batch-specific cell types with few cells, and a multi-batches dataset with imbalanced batches and batch-specific cell types. Moreover, it best preserves the relative gene expression between cell types, yielding superior differential expression analysis results. Finally, in a more complex setting of a Leukemia cohort, our method preserved most of the underlying biological information for each patient while aligning the batches, improving the clustering metrics in the aggregated dataset.

BR55 Ultrasound Molecular Imaging of Clear Cell Renal Cell Carcinoma Reflects Tumor Vascular Expression of VEGFR-2 in a Patient-Derived Xenograft Model.

Standard imaging cannot reliably predict the nature of renal tumors. Among malignant renal tumors, clear cell renal cell carcinoma (ccRCC) is the most common histological subtype, in which the vascular endothelial growth factor 2 (VEGFR-2) is highly expressed in the vascular endothelium. BR55, a contrast agent for ultrasound imaging, consists of gas-core lipid microbubbles that specifically target and bind to the extracellular portion of the VEGFR-2. The specific information provided by ultrasound molecular imaging (USMI) using BR55 was compared with the vascular tumor expression of the VEGFR-2 by immunohistochemical (IHC) staining in a preclinical model of ccRCC. Patients' ccRCCs were orthotopically grafted onto Nod-Scid-Gamma (NSG) mice to generate patient-derived xenografts (PdX). Mice were divided into four groups to receive either vehicle or axitinib an amount of 2, 7.5 or 15 mg/kg twice daily. Perfusion parameters and the BR55 ultrasound contrast signal on PdX renal tumors were analyzed at D0, D1, D3, D7 and D11, and compared with IHC staining for the VEGFR-2 and CD34. Significant Pearson correlation coefficients were observed between the area under the curve (AUC) and the CD34 (0.84, p < 10-4), and between the VEGFR-2-specific signal obtained by USMI and IHC (0.72, p < 10-4). USMI with BR55 could provide instant, quantitative information on tumor VEGFR-2 expression to characterize renal masses non-invasively.

Synergic prognostic value of 3D CT scan subcutaneous fat and muscle masses for immunotherapy-treated cancer.

BACKGROUND: Our aim was to explore the prognostic value of anthropometric parameters in a large population of patients treated with immunotherapy. METHODS: We retrospectively included 623 patients with advanced non-small cell lung cancer (NSCLC) (n=318) or melanoma (n=305) treated by an immune-checkpoint-inhibitor having a pretreatment (thorax-)abdomen-pelvis CT scan. An external validation cohort of 55 patients with NSCLC was used. Anthropometric parameters were measured three-dimensionally (3D) by a deep learning software (Anthropometer3DNet) allowing an automatic multislice measurement of lean body mass, fat body mass (FBM), muscle body mass (MBM), visceral fat mass (VFM) and sub-cutaneous fat mass (SFM). Body mass index (BMI) and weight loss (WL) were also retrieved. Receiver operator characteristic (ROC) curve analysis was performed and overall survival was calculated using Kaplan-Meier (KM) curve and Cox regression analysis. RESULTS: In the overall cohort, 1-year mortality rate was 0.496 (95% CI: 0.457 to 0.537) for 309 events and 5-year mortality rate was 0.196 (95% CI: 0.165 to 0.233) for 477 events. In the univariate Kaplan-Meier analysis, prognosis was worse (p<0.001) for patients with low SFM (<3.95 kg/m2), low FBM (<3.26 kg/m2), low VFM (<0.91 kg/m2), low MBM (<5.85 kg/m2) and low BMI (<24.97 kg/m2). The same parameters were significant in the Cox univariate analysis (p<0.001) and, in the multivariate stepwise Cox analysis, the significant parameters were MBM (p<0.0001), SFM (0.013) and WL (0.0003). In subanalyses according to the type of cancer, all body composition parameters were statistically significant for NSCLC in ROC, KM and Cox univariate analysis while, for melanoma, none of them, except MBM, was statistically significant. In multivariate Cox analysis, the significant parameters for NSCLC were MBM (HR=0.81, p=0.0002), SFM (HR=0.94, p=0.02) and WL (HR=1.06, p=0.004). For NSCLC, a KM analysis combining SFM and MBM was able to separate the population in three categories with the worse prognostic for the patients with both low SFM (<5.22 kg/m2) and MBM (<6.86 kg/m2) (p<0001). On the external validation cohort, combination of low SFM and low MBM was pejorative with 63% of mortality at 1 year versus 25% (p=0.0029). CONCLUSIONS: 3D measured low SFM and MBM are significant prognosis factors of NSCLC treated by immune checkpoint inhibitors and can be combined to improve the prognostic value.

A biology-driven deep generative model for cell-type annotation in cytometry.

Cytometry enables precise single-cell phenotyping within heterogeneous populations. These cell types are traditionally annotated via manual gating, but this method lacks reproducibility and sensitivity to batch effect. Also, the most recent cytometers-spectral flow or mass cytometers-create rich and high-dimensional data whose analysis via manual gating becomes challenging and time-consuming. To tackle these limitations, we introduce Scyan https://github.com/MICS-Lab/scyan, a Single-cell Cytometry Annotation Network that automatically annotates cell types using only prior expert knowledge about the cytometry panel. For this, it uses a normalizing flow-a type of deep generative model-that maps protein expressions into a biologically relevant latent space. We demonstrate that Scyan significantly outperforms the related state-of-the-art models on multiple public datasets while being faster and interpretable. In addition, Scyan overcomes several complementary tasks, such as batch-effect correction, debarcoding and population discovery. Overall, this model accelerates and eases cell population characterization, quantification and discovery in cytometry.

Optimizing IFN Alpha Therapy against Myeloproliferative Neoplasms.

Myeloproliferative neoplasms (MPNs) are hematologic malignancies that result from acquired driver mutations in hematopoietic stem cells (HSCs), causing overproduction of blood cells and an increased risk of thrombohemorrhagic events. The most common MPN driver mutation affects the JAK2 gene (JAK2V617F ). Interferon alpha (IFNα) is a promising treatment against MPNs by inducing a hematologic response and molecular remission for some patients. Mathematical models have been proposed to describe how IFNα targets mutated HSCs, indicating that a minimal dose is necessary for long-term remission. This study aims to determine a personalized treatment strategy. First, we show the capacity of an existing model to predict cell dynamics for new patients from data that can be easily obtained in clinic. Then, we study different treatment scenarios in silico for three patients, considering potential IFNα dose-toxicity relations. We assess when the treatment should be interrupted depending on the response, the patient's age, and the inferred development of the malignant clone without IFNα We find that an optimal strategy would be to treat patients with a constant dose so that treatment could be interrupted as quickly as possible. Higher doses result in earlier discontinuation but also higher toxicity. Without knowledge of the dose-toxicity relationship, trade-off strategies can be found for each patient. A compromise strategy is to treat patients with medium doses (60-120 µg/week) for 10-15 years. Altogether, this work demonstrates how a mathematical model calibrated from real data can help build a clinical decision-support tool to optimize long-term IFNα therapy for MPN patients. SIGNIFICANCE STATEMENT: Myeloproliferative neoplasms (MPNs) are chronic blood cancers. Interferon alpha (IFNα) is a promising treatment with the potential to induce a molecular response by targeting mutated hematopoietic stem cells. MPN patients are treated over several years, and there is a lack of knowledge concerning the posology strategy and the best timing for interrupting therapy. The study opens avenues for rationalizing how to treat MPN patients with IFNα over several years, promoting a more personalized approach to treatment.

New calibration setup for quantitative DCE-US imaging protocol: Toward standardization.

BACKGROUND: The DCE-US (Dynamic Contrast-Enhanced Ultrasonography) imaging protocol predicts the vascular modifications compared with Response Evaluation Criteria in Solid Tumors (RECIST) based mainly on morphological changes. A quantitative biomarker has been validated through the DCE-US multi-centric study for early monitoring of the efficiency of anti-angiogenic cancer treatments. In this context, the question of transposing the use of this biomarker to other types of ultrasound scanners, probes and settings has arisen to maintain the follow-up of patients under anti-angiogenic treatments. As a consequence, radiologists encounter standardization issues between the different generations of ultrasound scanners to perform quantitative imaging protocols. PURPOSE: The aim of this study was to develop a new calibration setup to transpose the DCE-US imaging protocol to the new generation of ultrasound scanners using both abdominal and linear probes. METHODS: This calibration method has been designed to be easily reproducible and optimized, reducing the time required and cost incurred. It is based on an original set-up that includes using a concentration splitter to measure the variation of the harmonic signal intensity, obtained from the Area Under the time-intensity Curve (AUC) as a function of various contrast-agent concentrations. The splitter provided four different concentrations simultaneously ranging from 12.5% to 100% of the initial concentration of the SonoVue contrast agent (Bracco Imaging S.p.A., Milan, Italy), therefore, measuring four AUCs in a single injection. The plot of the AUC as a function of the four contrast agent concentrations represents the intensity variation of the harmonic signal: the slope being the calibration parameter. The standardization through this method implied that both generations of ultrasound scanners had to have the same slopes to be considered as calibrated. This method was tested on two ultrasound scanners from the same manufacturer (Aplio500, Aplioi900, Canon Medical Systems, Tokyo, Japan). The Aplio500 used the settings defined by the initial multicenter DCE-US study. The Mechanical Index (MI) and the Color Gain (CG) of the Aplioi900 have been adjusted to match those of the Aplio500. The reliability of the new setup was evaluated in terms of measurement repeatability, and reproducibility with the agreement between the measurements obtained once the two ultrasound scanners were calibrated. RESULTS: The new setup provided excellent repeatability measurements with a value of 96.8%. Once the two ultrasound scanners have been calibrated for both types of probes, the reproducibility was excellent with the agreement between their respective quantitative measurement was at the lowest 95.4% and at the best 98.8%. The settings of the Aplioi900 (Canon Medical Systems) were adjusted to match those of the Aplio500 (Canon Medical Systems) and these validated settings were for the abdominal probe: MI = 0.13 and CG = 34 dB; and for the linear probe: MI = 0.10 and CG = 38 dB. CONCLUSION: This new calibration setup provided reliable measurements and enabled the rapid transfer and the use of the DCE-US imaging protocol on new ultrasound scanners, thus permitting a continuation of the therapeutic evaluation of patients through quantitative imaging.

Integrative Pan-Cancer Genomic and Transcriptomic Analyses of Refractory Metastatic Cancer.

Metastatic relapse after treatment is the leading cause of cancer mortality, and known resistance mechanisms are missing for most treatments administered to patients. To bridge this gap, we analyze a pan-cancer cohort (META-PRISM) of 1,031 refractory metastatic tumors profiled via whole-exome and transcriptome sequencing. META-PRISM tumors, particularly prostate, bladder, and pancreatic types, displayed the most transformed genomes compared with primary untreated tumors. Standard-of-care resistance biomarkers were identified only in lung and colon cancers-9.6% of META-PRISM tumors, indicating that too few resistance mechanisms have received clinical validation. In contrast, we verified the enrichment of multiple investigational and hypothetical resistance mechanisms in treated compared with nontreated patients, thereby confirming their putative role in treatment resistance. Additionally, we demonstrated that molecular markers improve 6-month survival prediction, particularly in patients with advanced breast cancer. Our analysis establishes the utility of the META-PRISM cohort for investigating resistance mechanisms and performing predictive analyses in cancer. SIGNIFICANCE: This study highlights the paucity of standard-of-care markers that explain treatment resistance and the promise of investigational and hypothetical markers awaiting further validation. It also demonstrates the utility of molecular profiling in advanced-stage cancers, particularly breast cancer, to improve the survival prediction and assess eligibility to phase I clinical trials. This article is highlighted in the In This Issue feature, p. 1027.

CustOmics: A versatile deep-learning based strategy for multi-omics integration.

The availability of patient cohorts with several types of omics data opens new perspectives for exploring the disease's underlying biological processes and developing predictive models. It also comes with new challenges in computational biology in terms of integrating high-dimensional and heterogeneous data in a fashion that captures the interrelationships between multiple genes and their functions. Deep learning methods offer promising perspectives for integrating multi-omics data. In this paper, we review the existing integration strategies based on autoencoders and propose a new customizable one whose principle relies on a two-phase approach. In the first phase, we adapt the training to each data source independently before learning cross-modality interactions in the second phase. By taking into account each source's singularity, we show that this approach succeeds at taking advantage of all the sources more efficiently than other strategies. Moreover, by adapting our architecture to the computation of Shapley additive explanations, our model can provide interpretable results in a multi-source setting. Using multiple omics sources from different TCGA cohorts, we demonstrate the performance of the proposed method for cancer on test cases for several tasks, such as the classification of tumor types and breast cancer subtypes, as well as survival outcome prediction. We show through our experiments the great performances of our architecture on seven different datasets with various sizes and provide some interpretations of the results obtained. Our code is available on (https://github.com/HakimBenkirane/CustOmics).

Better than RECIST and Faster than iRECIST: Defining the Immunotherapy Progression Decision Score to Better Manage Progressive Tumors on Immunotherapy.

PURPOSE: The objective of the study is to propose the immunotherapy progression decision (iPD) score, a practical tool based on patient features that are available at the first evaluation of immunotherapy treatment, to help oncologists decide whether to continue the treatment or switch rapidly to another therapeutic line when facing a progressive disease patient at the first evaluation. EXPERIMENTAL DESIGN: This retrospective study included 107 patients with progressive disease at first evaluation according to RECIST 1.1. Clinical, radiological, and biological data at baseline and first evaluation were analyzed. An external validation set consisting of 31 patients with similar baseline characteristics was used for the validation of the score. RESULTS: Variables were analyzed in a univariate study. The iPD score was constructed using only independent variables, each considered as a worsening factor for the survival of patients. The patients were stratified in three groups: good prognosis (GP), poor prognosis (PP), and critical prognosis (CP). Each group showed significantly different survivals (GP: 11.4, PP: 4.4, CP: 2.3 months median overall survival, P < 0.001, log-rank test). Moreover, the iPD score was able to detect the pseudoprogressors better than other scores. On the validation set, CP patients had significantly worse survival than PP and GP patients (P < 0.05, log-rank test). CONCLUSIONS: The iPD score provides oncologists with a new evaluation, computable at first progression, to decide whether treatment should be continued (for the GP group), or immediately changed for the PP and CP groups. Further validation on larger cohorts is needed to prove its efficacy in clinical practice.

Machine learning identifies a profile of inadequate responder to methotrexate in rheumatoid arthritis.

OBJECTIVES: Around 30% of patients with RA have an inadequate response to MTX. We aimed to use routine clinical and biological data to build machine learning models predicting EULAR inadequate response to MTX and to identify simple predictive biomarkers. METHODS: Models were trained on RA patients fulfilling the 2010 ACR/EULAR criteria from the ESPOIR and Leiden EAC cohorts to predict the EULAR response at 9 months (± 6 months). Several models were compared on the training set using the AUROC. The best model was evaluated on an external validation cohort (tREACH). The model's predictions were explained using Shapley values to extract a biomarker of inadequate response. RESULTS: We included 493 therapeutic sequences from ESPOIR, 239 from EAC and 138 from tREACH. The model selected DAS28, Lymphocytes, Creatininemia, Leucocytes, AST, ALT, swollen joint count and corticosteroid co-treatment as predictors. The model reached an AUROC of 0.72 [95% CI (0.63, 0.80)] on the external validation set, where 70% of patients were responders to MTX. Patients predicted as inadequate responders had only 38% [95% CI (20%, 58%)] chance to respond and using the algorithm to decide to initiate MTX would decrease inadequate-response rate from 30% to 23% [95% CI: (17%, 29%)]. A biomarker was identified in patients with moderate or high activity (DAS28 > 3.2): patients with a lymphocyte count superior to 2000 cells/mm3 are significantly less likely to respond. CONCLUSION: Our study highlights the usefulness of machine learning in unveiling subgroups of inadequate responders to MTX to guide new therapeutic strategies. Further work is needed to validate this approach.

Predicting immunotherapy outcomes in patients with MSI tumors using NLR and CT global tumor volume.

Background: Anti-PD-(L)1 treatment is indicated for patients with mismatch repair-deficient (MMRD) tumors, regardless of tumor origin. However, the response rate is highly heterogeneous across MMRD tumors. The objective of the study is to find a score that predicts anti-PD-(L)1 response in patients with MMRD tumors. Methods: Sixty-one patients with various origin of MMRD tumors and treated with anti-PD-(L)1 were retrospectively included in this study. An expert radiologist annotated all tumors present at the baseline and first evaluation CT-scans for all the patients by circumscribing them on their largest axial axis (single slice), allowing us to compute an approximation of their tumor volume. In total, 2120 lesions were annotated, which led to the computation of the total tumor volume for each patient. The RECIST sum of target lesions' diameters and neutrophile-to-lymphocyte (NLR) were also reported at both examinations. These parameters were determined at baseline and first evaluation and the variation between the first evaluation and baseline was calculated, to determine a comprehensive score for overall survival (OS) and progression-free survival (PFS). Results: Total tumor volume at baseline was found to be significantly correlated to the OS (p-value: 0.005) and to the PFS (p-value:<0.001). The variation of the RECIST sum of target lesions' diameters, total tumor volume and NLR were found to be significantly associated to the OS (p-values:<0.001, 0.006,<0.001 respectively) and to the PFS (<0.001,<0.001, 0.007 respectively). The concordance score combining total tumor volume and NLR variation was better at stratifying patients compared to the tumor volume or NLR taken individually according to the OS (pairwise log-rank test p-values: 0.033,<0.001, 0.002) and PFS (pairwise log-rank test p-values: 0.041,<0.001, 0.003). Conclusion: Total tumor volume appears to be a prognostic biomarker of anti-PD-(L)1 response to immunotherapy in metastatic patients with MMRD tumors. Combining tumor volume and NLR with a simple concordance score stratifies patients well according to their survival and offers a good predictive measure of response to immunotherapy.

Inferring the initiation and development of myeloproliferative neoplasms.

The developmental history of blood cancer begins with mutation acquisition and the resulting malignant clone expansion. The two most prevalent driver mutations found in myeloproliferative neoplasms-JAK2V617F and CALRm-occur in hematopoietic stem cells, which are highly complex to observe in vivo. To circumvent this difficulty, we propose a method relying on mathematical modeling and statistical inference to determine disease initiation and dynamics. Our findings suggest that CALRm mutations tend to occur later in life than JAK2V617F. Our results confirm the higher proliferative advantage of the CALRm malignant clone compared to JAK2V617F. Furthermore, we illustrate how mathematical modeling and Bayesian inference can be used for setting up early screening strategies.

Machine learning predicts response to TNF inhibitors in rheumatoid arthritis: results on the ESPOIR and ABIRISK cohorts.

OBJECTIVES: Around 30% of patients with rheumatoid arthritis (RA) do not respond to tumour necrosis factor inhibitors (TNFi). We aimed to predict patient response to TNFi using machine learning on simple clinical and biological data. METHODS: We used data from the RA ESPOIR cohort to train our models. The endpoints were the EULAR response and the change in Disease Activity Score (DAS28). We compared the performances of multiple models (linear regression, random forest, XGBoost and CatBoost) on the training set and cross-validated them using the area under the receiver operating characteristic curve (AUROC) or the mean squared error. The best model was then evaluated on a replication cohort (ABIRISK). RESULTS: We included 161 patients from ESPOIR and 118 patients from ABIRISK. The key selected features were DAS28, lymphocytes, ALT (aspartate aminotransferase), neutrophils, age, weight, and smoking status. When predicting EULAR response, CatBoost achieved the best performances of the four tested models. It reached an AUROC of 0.72 (0.68-0.73) on the train set (ESPOIR). Better results were obtained on the train set when etanercept and monoclonal antibodies were analysed separately. On the test set (ABIRISK), these models respectively achieved on AUROC of 0.70 (0.57-0.82) and 0.71 (0.55-0.86). Two decision thresholds were tested. The first prioritised a high confidence in identifying responders and yielded a confidence up to 90% for predicting response. The second prioritised a high confidence in identifying inadequate responders and yielded a confidence up to 70% for predicting non-response. The change in DAS28 was predicted with an average error of 1.1 DAS28 points. CONCLUSION: The machine learning models developed allowed predicting patient response to TNFi exclusively using data available in clinical routine.

Prediction of Early Response to Immunotherapy: DCE-US as a New Biomarker.

PURPOSE: The objective of our study is to propose fast, cost-effective, convenient, and effective biomarkers using the perfusion parameters from dynamic contrast-enhanced ultrasound (DCE-US) for the evaluation of immune checkpoint inhibitors (ICI) early response. METHODS: The retrospective cohort used in this study included 63 patients with metastatic cancer eligible for immunotherapy. DCE-US was performed at baseline, day 8 (D8), and day 21 (D21) after treatment onset. A tumor perfusion curve was modeled on these three dates, and change in the seven perfusion parameters was measured between baseline, D8, and D21. These perfusion parameters were studied to show the impact of their variation on the overall survival (OS). RESULTS: After the removal of missing or suboptimal DCE-US, the Baseline-D8, the Baseline-D21, and the D8-D21 groups included 37, 53, and 33 patients, respectively. A decrease of more than 45% in the area under the perfusion curve (AUC) between baseline and D21 was significantly associated with better OS (p = 0.0114). A decrease of any amount in the AUC between D8 and D21 was also significantly associated with better OS (p = 0.0370). CONCLUSION: AUC from DCE-US looks to be a promising new biomarker for fast, effective, and convenient immunotherapy response evaluation.

Representation Learning for the Clustering of Multi-Omics Data.

The integration of several sources of data for the identification of subtypes of diseases has gained attention over the past few years. The heterogeneity and the high dimensions of the data sets calls for an adequate representation of the data. We summarize the field of representation learning for the multi-omics clustering problem and we investigate several techniques to learn relevant combined representations, using methods from group factor analysis (PCA, MFA and extensions) and from machine learning with autoencoders. We highlight the importance of appropriately designing and training the latter, notably with a novel combination of a disjointed deep autoencoder (DDAE) architecture and a layer-wise reconstruction loss. These different representations can then be clustered to identify biologically meaningful clusters of patients. We provide a unifying framework for model comparison between statistical and deep learning approaches with the introduction of a new weighted internal clustering index that evaluates how well the clustering information is retained from each source, favoring contributions from all data sets. We apply our methodology to two case studies for which previous works of integrative clustering exist, TCGA Breast Cancer and TARGET Neuroblastoma, and show how our method can yield good and well-balanced clusters across the different data sources.

Inferring the dynamics of mutated hematopoietic stem and progenitor cells induced by IFNα in myeloproliferative neoplasms.

Classical BCR-ABL-negative myeloproliferative neoplasms (MPNs) are clonal disorders of hematopoietic stem cells (HSCs) caused mainly by recurrent mutations in genes encoding JAK2 (JAK2), calreticulin (CALR), or the thrombopoietin receptor (MPL). Interferon α (IFNα) has demonstrated some efficacy in inducing molecular remission in MPNs. To determine factors that influence molecular response rate, we evaluated the long-term molecular efficacy of IFNα in patients with MPN by monitoring the fate of cells carrying driver mutations in a prospective observational and longitudinal study of 48 patients over more than 5 years. We measured the clonal architecture of early and late hematopoietic progenitors (84 845 measurements) and the global variant allele frequency in mature cells (409 measurements) several times per year. Using mathematical modeling and hierarchical Bayesian inference, we further inferred the dynamics of IFNα-targeted mutated HSCs. Our data support the hypothesis that IFNα targets JAK2V617F HSCs by inducing their exit from quiescence and differentiation into progenitors. Our observations indicate that treatment efficacy is higher in homozygous than heterozygous JAK2V617F HSCs and increases with high IFNα dose in heterozygous JAK2V617F HSCs. We also found that the molecular responses of CALRm HSCs to IFNα were heterogeneous, varying between type 1 and type 2 CALRm, and a high dose of IFNα correlates with worse outcomes. Our work indicates that the long-term molecular efficacy of IFNα implies an HSC exhaustion mechanism and depends on both the driver mutation type and IFNα dose.

Ultrasound Molecular Imaging of Renal Cell Carcinoma: VEGFR targeted therapy monitored with VEGFR1 and FSHR targeted microbubbles.

Recent treatment developments for metastatic renal cell carcinoma offer combinations of immunotherapies or immunotherapy associated with tyrosine kinase inhibitors (TKI). There is currently no argument to choose one solution or another. Easy-to-use markers to assess longitudinal responses to TKI are necessary to determine when to switch to immunotherapies. These new markers will enable an earlier adaptation of therapeutic strategy in order to prevent tumor development, unnecessary toxicity and financial costs. This study evaluates the potential of ultrasound molecular imaging to track the response to sunitinib in a clear cell renal carcinoma model (ccRCC). We used a patient-derived xenograft model for this imaging study. Mice harboring human ccRCC were randomized for sunitinib treatment vs. control. The tumors were imaged at days 0, 7, 14 and 28 with ultrasound molecular imaging. Signal enhancement was quantified and compared between the two groups after injections of non-targeted microbubbles and microbubbles targeting VEGFR1 and FSHR. The tumor growth of the sunitinib group was significantly slower. There was a significantly lower expression of both VEGFR-1 and FSHR molecular ultrasound imaging signals in the sunitinib group at all times of treatment (Days 7, 14 and 28). These results confirm the study hypothesis. There was no significant difference between the 2 groups for the non-targeted microbubble ultrasound signal. This study demonstrated for the first time the potential of VEGFR1 and FSHR, by ultrasound-based molecular imaging, to follow-up the longitudinal response to sunitinib in ccRCC. These results should trigger developments for clinical applications.

Association of metastatic pattern and molecular status in stage IV non-small cell lung cancer adenocarcinoma.

OBJECTIVES: The aim of our study was to investigate the association between driver oncogene alterations and metastatic patterns on imaging assessment, in a large cohort of metastatic lung adenocarcinoma patients. METHODS: From January 2010 to May 2017, 550 patients with stage IV lung adenocarcinoma with molecular analysis were studied retrospectively including 135 EGFR-mutated, 81 ALK-rearrangement, 47 BRAF-mutated, 141 KRAS-mutated, and 146 negative tumors for these 4 mutations (4N). After review of the complete imaging report by two radiologists (junior and senior) to identify metastatic sites, univariate correlation analyzes were performed. RESULTS: We found differences in metastatic tropism depending on the molecular alteration type when compared with the non-mutated 4N group: in the EGFR group, pleural metastases were more frequent (32% versus 20%; p = 0.021), and adrenal and node metastases less common (6% versus 23%; p < 0.001 and 11% versus 23%; p = 0.011). In the ALK group, there were more brain and lung metastases (respectively 42% versus 29%; p = 0.043 and 37% versus 24%; p = 0.037). In the BRAF group, pleural and pericardial metastases were more common (respectively 47% versus 20%; p < 0.001 and 11% versus 3%; p = 0.04) and bone metastases were rarer (21% versus 42%; p = 0.011). Lymphangitis was more frequent in EGFR, ALK, and BRAF groups (respectively 6%, 7%, and 15% versus 1%); p = 0.016; p = 0.009; and p < 0.001. CONCLUSION: The application of these correlations between molecular status and metastatic tropism in clinical practice may lead to earlier and more accurate identification of patients for targeted therapy. KEY POINTS: • Bone and brain metastasis are the most common organs involved in lung adenocarcinoma but the relative incidence of each metastatic site depends on the molecular alteration. • EGFR-mutated tumors preferentially spread to the pleura and less commonly to adrenals, ALK-rearrangement tumors usually spread to the brain and the lungs, whereas BRAF-mutated tumors are unlikely to spread to bones and have a serous (pericardial ad pleural) tropism. • These correlations could help in the clinical management of patients with metastatic lung adenocarcinoma.